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Abstract 

A comparative measurement of timing properties of magnetic-resistant 
fine mesli R7761-70 and ordinary fast R2083 photomultipliers is pre- 
sented togetlier with preliminary results on the operation of R7761-70 
PMs in magnetic field up to 1100 Gauss. The results were obtained 
using the proton beam of the MC50 Cyclotron of Korea Institute of 
Radilogical and Medical Sciences. 

The ratio of the effective R7761-70 and R2083 TOF (or timing) 
resolutions was extracted by using two different methods. The results 
are 1.05 ± 0.066 and 1.07 ± 0.062. The gain of R7761-70 PMs is 
not affected by magnetic field. The R7761-70 TOF/timing resolution 
becomes ~ 8% better at 1100 Gauss if the external field is oriented 
parallel to the PM axis. The results prove the advantages of the 
design of the CLAS12 Central Time-of- flight system with fine-mesh 
photomultipliers in comparison with the "conservative" design based 
on ordinary R2083 PMs and long bent light guides. 

1 Introduction 

Initially the Central Time-of-Flight system (CTOF) was considered as a bar- 
rel made of ~ 50 scintillator bars viewed by fast Hamamatsu R2083 pho- 
tomultipliers (PMs) through ~ 1.5 — 1.7 m long bent hght-guides[Il |2] (left 
panel of Fig. [1]). The light guides are to deliver scintillation light to the 
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regions outside the central solenoid where the magnetic field drops down to 
300 - 1000 Gauss. Being properly shielded, ordinary photomultipliers are 
expected to operate in this field. 

Problems of this design are obvious: 

- The long and bent light guides would deliver a small portion of scintilla- 
tion light to PMs thus deteriorating the time-of-fiight (TOF) resolution; 

- The mechanic construction would be rather complicate, ugly, and frag- 
ile. It would interfere with other CLAS12 sub-detectors; 

- ~ 100 long light guides and ~ 100 magnetic shields would significantly 
increase the overall CTOF weight and costs; 

The Nuclear Physics Group of Kyungpook National University, in col- 
laboration with Inisitute for Nuclear Research, Moscow, suggests another 
solution based on magnetic-resistant fine-mesh photomultipliers. Fine-mesh 
photomultipliers can operate in magnetic field up to 1.5 Tesla. They could 
be placed closer to scintillator bars at positions where the magnetic field 
is 0.3 - 1 Tesla (right panel of Fig. [T]). The light guides would be much 
shorter, ~ 0.50 — 0.8 m long, and not bent. No shields would be needed. 
The CTOF assembly would be simpler, more reliable, and less expensive. A 
critical question is whether the acceptable TOF resolution could be achieved 
with fine-mesh of photomultipliers in comparison with R2083 PMs, and how 
it would be affected by the magnetic field 

There was some skepticism regarding fine-mesh photomultipliers. The 
authors of Ref. [2] wrote that "...While the transition time spread (TTS) 
is compatible to TTS of R2083, its anode rise time [Rem. of a fine-mesh 
photomultiplier) is 4 times larger. Therefore it is unlikely it may be used 
for precise timing...". Our consideration is different: the real PM anode rise 
time in a scintillation counter depends also on the light decay constant of a 
scintillator material. For Bicron-408 it is 2.5 ns - is much longer than the 
internal R2083 rise time. The most critical parameter is the transit time 
spread (TTS) which is in fact the variation of electron propagation time 
inside a photomultiplier. The TTS of fine-mesh R7761-70 photomultipliers 
is 350 ps. It is better than that of R2083 PMs (370 ps) g]. 

The operation of fine-mesh photomultipliers in magnetic field was inves- 
tigated in Ref. |5l|6l[7]. The outcome of those studies was encouraging: the 
timing performance of fine-mesh photomultipliers is good and is not affected 
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Figure 1: Two versions of the CTOF design. On the left: "conservative" 
design with long bent light guides and ordinary R2083 photomultipliers. On 
the right: CTOF design with fine- mesh photomultipliers with short straight 
light guides. 

by magnetic field up to ~ 0.8 Tesla. The authors employed fast laser light 
pulses as a tool to study PM operation. This technique is quite convenient 
but it doesn't allow to estimate the CTOF resolution and to judge between 
two CTOF designs. 

Until recently the cosmic-ray tracking was believed to be a main tool for 
the CTOF R&Ds [21 E]- The method uses three stacked parallel equidistant 
counters viewed by six phototubes under study. If a cosmic-ray muon crosses 
all three counters, the time (or coordinate) of a scintillation in the middle 
counter is equal to a half of the sum of times/coordinates in the top and 
bottom counter^. This leads to the following relation 

T = —{tmidl + tmid2) ~ ^ (^topl + ttop2 + tbotl + ^60*2) = iX) 

where ttopi----tbot2 denote 6 PM signal arriving times derived from TDC read- 
outs. From Eq. [T]one may deduce the effective timing resolution in each PM 
channel 

2 

(TpM = yl^(^) (2) 

In practice the PM timing resolution is extracted from the width of the peak 
in the spectrum of r. 

By employing this method, the KNU group measured the effective R7761- 
70 timing resolution cpMmYei-rofi ~ 52 ps and found it similar to that of 

^Detailed description of this method is available in [3] 
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R2083 <7pM2083fi ~ 53 ps [S]. The second result is in good agreement with 
the previous measurement [5] (TpM2083fi = 52ps ± 1 ps. 

A disadvantage of this technique is that the measured PM arriving times 
ttopi----tbot2 are spoiled by time walks of constant-fraction discriminators (CFD) 
Time walks depend on pulse shapes and especially on pulse heights. The CFD 
Phillips 715 and Ortec935 used in Ref . [21 El [H] , generate time walks up to 
±50ps . Due to light attenuation the pulse heights of muons events depend 
on hit coordinates and track angles. The variation of the corresponding time 
walks is comparable with the measured PM timing resolution. That is why 
the results obtained with cosmic rays, require sophisticated corrections and 
are cut-, electronic-, and analysis-dependentH. They can be used only for 
rough preliminary estimates of the expected TOF resolution. 

In this Note, we present another method to study the operation of scin- 
tillation counters and photomultipliers using a well-coUimated proton beam. 
The method has been implemented at the MC50 Cyclotron of Korea Institute 
of Radiological and Medical Sciences. Simultaneously we report the measure- 
ment of the TOF resolution of a plastic-scintillation counter equipped with 
fine-mesh Hamamatsu R7761-70 photomultipliers and compare it with that 
obtained with ordinary fast Hamamatsu R2083 PMs. We also report our 
first results on the operation of fine-mesh R7761-70 photomultipliers in the 
magnetic field up to 1100 Gauss. 

2 Fine-mesh photomultipliers 

Fine-mesh photomultipliers have been developed for high magnetic-field ap- 
plications. Their dynode system has a structure of fine-mesh electrodes 
stacked in close proximity. Such dynodes provide an improved pulse linearity 
and resistance to external magnetic field [Ij. 

The survey of properties of fine- mesh photomultipliers is given in Table [1] 
together with those of R2083 PMs. In general, the timing characteristics of 
the fine- mesh photomultipliers are worse: the anode rise time varies from 2.1 
to 2.7 ns. For R2083 PMs this number is 0.7 ns. However the R7761-70 and 
R5505-70 transit time spreads are better than that of R2083 PMs. 

Due to geometrical dimensions only R7761-70 and R5924-70 PMs are 
suitable for CTOF. Among them, R7761-70 PMs were chosen for initial tests 
because 

^The influence of time walks on cosmic-ray results is discussed in detail in Ref. [5]. 
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Phototube 


R2083 


R5505-70 


R7761-70 


R5924-70 


Dynode system 


ordinary 


fine-mesh 


fine-mesh 


fine-mesh 


Photocathode dia (mm) 


39 


17.5 


27 


39 


Photocathode type 


bialkali 


bialkali 


bialkali 


bialkali 


Anode sensitivity (A/lm) 


200 


40 


800 


700 


Anode rise time (ns) 


0.7 


2.1 


2.5 


2.5 


Transit time (ns) 


16 


5.6 


7.5 


9.5 


Transit time 


0.37 


0.35 


0.35 


0.44 



spread (ns) 



Table 1: Characteristics of fine- mesh photomultipliers in comparison with 
R2083. 

i) Their timing characteristics are better than those of R5924-70; 

ii) Small dimensions allow us to use them at the CTOF upstream ends; 

iii) These phototubes are essentially cheaper than R5924-70 PMs. 
Hamamatsu Photonics offers H8409-70 assemblies. A H8409-70 assembly 

consists of a R7761-70 PM, a voltage divider designed for positive HV supply, 
and a phototube housing. A positive divider necessarily includes a capacitor 
in the anode circuit that separates the voltage supply from the anode output. 
This capacitor may deteriorate timing properties and generate some level at 
high count rates. The geometric dimensions of H8409-70 assembly are larger 
than those of a single R7761-70 PM because of a phototube housing. 

We developed our own voltage divider designed for negative HV. The 
distribution of potentials in the dynode system was optimized to achieve the 
best rise time. R7761-70 PMs were optically attached to scintillator bars 
and wrapped round with isolation tape without any housing. This made it 
possible to minimize the counter dimensions and to reduce bores of solenoids 
used in magnetic-field measurements (Section [TTi). 

Typical R7761-70 and R2083 signals are shown in Fig. [21 They correspond 
to the detection of cosmic-ray muons in a 3 cm thick scintillator counter. The 
signals were obtained with HV ~ 22mV for R7761-70 PMs and HV ~ 2500 
V for R2083 PMs. The R7761-70 rise time ~ 3.2 ns. It is only slightly 
worse than that obtained with R2083PMs (~ 2.5 ns). The R7761-70 pulse 
height reaches 12V/50Ohm. Such high pulse heights assure the operation of 
fine-mesh PMs in magnetic field in which the PM gain might be lower. 
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Figure 2: Typical anode pulses from cosmic-ray muons obtained with R2083 
(left) and R7761-70(middle and right). 

3 Method and experimental setup 

If a long scintillator bar is viewed by two photomultipliers, the PM arriving 
times ti and ^2 corresponding to a particle-induced scintillation are defined 
by the following relations: 

ti = TOF + x/v + Canst = TOF + {L~x)/v + Const (3) 

where TOF is the time-of- flight of a particle from a certain point (target), 
X is the coordinate of a scintillation along the counter axis, L is the total 
length of a bar, v is the efficient speed of light propagation inside a bar, 
constants originate from cable and electronic delays. Therefore TOF and the 
coordinate can be derived from the PM times 

TOF = {ti+t2)/2 + Const (4) 

x = v{ti- t2) /2 + Const (5) 

The TOF resolution is 

cttof = (r{{ti + t2)/2) = a\ti) + a\t2) (6) 

where cr(ti) and (T{t2) are the effective timing resolutions in each PM channel. 
The variation of the time difference is equal to the TOF resolution 

a((ti - t2)/2) = ya^{t,)+a^{t2) = aroF (7) 

If a scintillation counter is irradiated by a narrow proton beam, the beam 
generates a peak in the distribution of events over the coordinate x. As 
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Figure 3: Experimental setup. 



it follows from Eq. [5], the coordinate distribution is equivalent to the time- 
difference spectrum {ti — t2)/2 scaled by factor v. The width of this peak in 
the time-difference spectrum is defined by the size of the beam spot Ax and 
by the timing performance of the counter 

(Tpeak = cr((ti - t2)/2 + Ax/v) = ^ a^oF + ^^i^^ / v)\ (8) 

where a{Ax/v) denotes rms in the coordinate distribution of events due to 
the finite beam dimension. For a point-like beam Ax ~ 

(7peak ~ CTtoF (9) 

A well-collimated proton beam allows to derive the TOF resolution by mea- 
suring the time-difference spectrum (ti — t2)/2. 

In the practical implementation of this method a counter made of the 
50 X 3 X 2 cm? Bicron-408 plastic scintillator bar was used. The bar was 
viewed by two photomultipliers under study. 

The counter was irradiated by the proton beam of the MC50 Cyclotron of 
Korea Institute of Radiological and Medical Sciences (KIRAMS). The beam 
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was collimated by a collimator made of 3 stacked together 7 mm thick steel 
plates. The plates had successively reducing holes of 5, 3, and 1 mm diameter 
in the first, second, and third plates respectively. The plate with the 1 mm 
dia hole was attached directly to the 2-cm side of the counter. The whole 
collimation assembly and the counter were fixed at the end of the beam pipe. 
As it was deduced from the Geant4 simulations, such collimator reduces the 
contamination of protons scattered on the collimator walls and minimizes 
the size of the beam spot. 

The PM times ti and ^2 and their pulse heights Ai and A2 were digitized 
by LeCroy 2228B TDCs and LeCroy2229B QDCs and were recorded on- 
line (Fig. [3]). One photomultiplier generated the common TDC START and 
triggered the acquisition. Both photomultipliers generated the STOP signals 
for TDCs. 

Typically, the beam current was set to 0.15 nA. At this current the count 
rate in the counter was ~ 2 x 10^ — 10^ Hz (Sect. [6]). The high voltages for 
fine-mesh R7761-70 PMs were set to low values of 1300 — 1350 V, in order 
to fit QDC ranges. 

4 Beam of the MC50 Cyclotron 

The MC50 Cyclotron of KIRAMS has been built in 1985 for medical, nuclear- 
physics, and biological applications. It can produce 20-50 MeV protons and 
deuterons. A neutron beam line is under construction. Beyond of a medical 
facility, there are three experimental hatches available for scientific research. 
The beam spectrum and profile in each hatch depend on the accelerator 
adjustment, the internal collimators, and other factors. They have to be 
adjusted for each experiment. 

The data reported here were taken as a sequence of short (~ 1 minute) 
measurements. In total there were 6 beam runs. Among them the data from 
three last runs are used in this Note. Each beam run lasted two-three days. It 
included apparatus installation at the MC50 beam line, calibrations, beam 
adjustment, and main measurements. The beam adjustment was a daily 
starting point and the most complicated and time-consuming procedure. The 
beam focusing and the position of the beam spot at the front end of the 
beam pipe were varied by the accelerator operator. The measured beam 
light-output spectrum was used as the criterion for the optimization of the 
beam quality. An example of this procedure is shown in Fig. HJ 8 measured 
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beam spectra correspond to 8 successive steps in the beam tuning. At the 
final point the set of beam parameters corresponding to the beam spectrum 
6 was chosen for the further running. 

The adjusted beam spectrum (the panels 5 and 6 of Fig. Hj) exhibits a 
main peak that corresponds to direct protons, and a tail at lower energies 
that is produced by particles which were either not properly accelerated or 
scattered inside the collimator. The position of the peak (i.e. the proton 
energy) depends on the operating parameters. Usually it was the same while 
in some cases the peak position was rather different. Such spectrum is shown 
in the panel 6 of Fig. HI There the peak position is nearly twice lower than 
its normal value. Other beam spectra from the runs used in this report 
are shown in Fig. [71 [HI The calibration of these spectra is explained in the 
section [51 



5 Birks' effect and light-output calibration 



Plastic scintillators do not respond linearly to the ionization density. Very 
dense ionization columns emit less light than that expected on the basis of 
^ for minimum ionizing particles. This non-linear response of scintillators 
is called Birks' effect [lOj. Due to Birks's effect, protons which stop inside 
a detector, produce less light per unit of deposited energy than relativistic 
minimum-ionizing particles. 

The semiempirical Birks' law is 



dL 

/ ax 

dE 



where L is the scintillator light production, Lq is the specific light production 
at low ionization densities (i.e. the light produced by a relativistic minimum- 
ionizing particle per a unit of deposited energy), x is a coordinate along the 
particle track inside a scintillator volume, and kb is Birks' constant which 
must be determined for each scintillator experimentally. An interpretation 
of Birks' effect was proposed by C.Chou pjj. He corrected Birks' formula as 

dT ^ 
where k}y and are adjustable constant. 
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Figure 4: Procedure of beam adjustment. 8 beam spectra correspond to 8 
successive steps in beam tuning. Calibration of the light output L is explained 
in the Section [5l 
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One can write Birks' law in the general form 



dL 

dx 



kp{E)Lo 



dE 
dx 



(12) 



where kp{E) is the Birks's coefficient that depends on E, on scintillator ma- 
terial, and on the type of a particle. 

Low-energy particles generate less light than relativistic minimum-ionizing 
(MI) particles. With the increase of the particle energy (/3 — >■ 1), the ratio of 
light production to deposited energy asymptotically approaches a constant 
value Lq which is the same for all particles. The latter is clear, for exam- 
ple, from the bi- dimensional plots light-output vs TOF for charged pions 
and protons obtained with the forward lead-scintillator TOF wall p3j at the 
GRAAL facility. It is convinient to assign this value to Lq = 1. In this case 
the energy deposited by a minimum-ionizing particle in a detector volume 
can be used as a measure for light output. 

TOF resolution depends on the number of photoelectrons produced at 
PM photocathodes. The latter is proportional to light output. The CTOF 
R&D requires to extrapolate the results obtained with the MC50 protons to 
those expected for fast minimum-ionizing particles. That is why the proton 
light output was calibrated by assigning it to the light output produced by 
high-energy cosmic-ray muons. 

The muon spectrum was measured just before and/or just after beam 
measurements. Schematic view of this measurement is shown in Fig. O The 
counter was disattached from the beam pipe and turned at 90° such that 
the effective counter thickness of 3 cm was the same for the muons and the 
beam protons. The count rate of muon events is essentially lower than the 
background in the experimental hutch. To reject this background, the second 
similar counter was placed 35 cm below the main counter. The coincidence of 
four signals from two counters was requested to trigger the acquisition. This 
made it possible to select on-line mostly those events in which a cosmic-ray 
muon passes both counters (Fig. [5]). 

Beam protons stop inside the counter and totally deposit their energies 
in the counter volume. Fast muons pass through the counter and act as 
minimum-ionizing particles. Their energy depositions depend on the track 
length inside the counter, i.e. on the angle between the muon trajectory and 
the counter surface. 

Two samples of the muon spectra are shown in Fig. [61 They contain 
the asymmetric peak with the maximum at ~ 6 MeV. This maximum corre- 
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Figure 5: Schematic two-side view of the cahbration measurement with 
cosmic-ray muons. 



sponds to the energy deposited by a minimum-ionizing particle crossing the 
counter in the perpedicular direction. The higher-energy part is formed by 
those muons which cross the counter at smaller angles. 

The light output can be derived from the QDCs readouts Ai and A2 after 
the pedestal subtraction as 

L = C^{Ai-pedi){A2-ped2) (13) 

where C is the calibration coefficients that relates QDC channels to the light 
output. The pedestals were measured by delaying the QDC gates for 200 ns. 
The calibration coefficient C was obtained by comparing the measured muon 
spectrum with the simulated spectrum of the energy deposited by muons 

(Fig. ED. 

As a cross-check, this calibration was applied to the data collected in 
the same beam run (January 6, 2009) with two counters. One counter was 
equipped with fine-mesh photomultipliers and the other with ordinary R2083 
PMs. The data were collected one-by-one keeping fixed the beam tuning. The 
muon spectrum with R7761-70 photomultipliers was collected just before the 
beam run while the one with R2083 PMs was taken immediately after. Both 
muon spectra are shown in Fig. [61 The calibration coefficients extracted from 
the muon spectra were then used to reconstruct the beam spectra. 
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Figure 6: Calibrated muon light-output spectra (Data from January 6, 2009). 
On the left: the spectrum obtained with fine-mesh R7761-70 photomultipli- 
ers. On the right: the spectrum obtained with ordinary R2083 photomulti- 
phers. Solid lines are simulated spectra of energy deposited by muons. 
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Figure 7: Beam spectra measured with the two different counters (Data from 
January 6 2009). On the left: the spectrum obtained with fine-mesh R7761- 
70 photomultipliers. On the right: the spectrum obtained with ordinary 
R2083 photomultipliers. 
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Figure 8: Beam spectra measured in October and November 2008. On the 
left: the spectrum obtained with fine-mesh R7761-70 photomuhiphers (Oc- 
tober 2008). On the right: the spectrum obtained with ordinary R2083 
photomuhiphers (November 2008). 

The resuhs are shown in Fig. [71 The beam peaks are located at the 
similar positions near 31.5 MeV. This proves the quality of the calibration 
procedure. The value of 31.5 MeV corresponds to the light output produced 
by a minimum-ionizing particle that deposit AE = 31.5 MeV in the counter 
volume. The beam protons stop inside the counter and deposit more energy, 
but, due to Birks' effect, generate the same light output. 

The position of the peak in the spectrum taken next day (January 7 
2009) is essentially lower, about ~ 17.5 MeV (the panel 6 of Fig. H]). This 
fact illustrates the importance and the infiuence of the MC50 beam tuning. 

The other beam spectra from the beam runs used in this report (October 
and November 2008) are shown in Fig. [HI The peaks are located at L ^ 
31.5 MeV. However the shapes of these spectra are different: the peak is 
wider and is smoothly continued by the low-energy tail. The reason for this 
difference still has to be understand. Our assumption is the internal MC50 
beam collimator was removed during October/November 2008 beam runs. 
The data from October/November 2008 were found suitable to retrieve the 
dependence of the time-of-flight resolution on the light output (Sect. [H]). 
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Figure 9: Beam spectra in logarithmic 
Fig. El 
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scale. The legend is the same as in 



The count rate in the counter depends on the beam current and tuning. In 
the measurements presented in this Note, the beam current was set to its 
lowest possible value 0.15 nA. The count rate can be estimated from the 
contamination of pileups in the beam spectrum. The probability that two 
pulses are encoded by QDC as one is 

Npileup _ 

Ppileup 77 wNqj. (14) 

where NpUeup is the number of pileup events in a beam spectrum, Ntot is the 
total number of events, w is the width of the QDC gates (100 ns in our case), 
and is the count rate in the counter. 

Fig. [9] shows a typical beam spectrum in the logarithmic scale. It contains 
events whose pulse heights are higher than the peak position. These are 
pileups. The fraction of such events typically varies from 0.2 to 1%. Following 
Eq. [131 one may estimate the count rate as 

A/'c^ ^ 2 X 10^ - l^fHz (15) 

At this count rate and the high voltage ~ 1300V the average anode 
current of R7761-70 photomultipliers was in the range of 1.6 — SyuA. Some 
tentative measurements were carried out at higher count rate/average anode 
current up to 5 x 10^ IIz/~ AQ^A. 
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It is worth noting that Ref. [5] quotes the maximum average anode current 
(typically lOOfiA for R5024-70 PMs ) which must not be exceeded. Beyond 
this limit the PM gain is expected to sharply drop down. For those two 
R7761-70 samples used in our tests, no reduction of the gain at the anode 
current up to was observed. More comprehensive tests at high count 

rates are planned for 2009. 



7 Getting TOF resolution 

An example of the measured [ti —t^)/^ spectrum is shown in Fig. [TOl The 
spectrum contains a narrow peak generated by the beam protons. The width 
of the peak can be extracted by means a Gaussian fit of this peak. However, 
because of the limited TDC resolution of ~ 47 ps/ch, the result is quite 
sensitive to the histogram binning. We found critical to set the width of the 
histogram bins equal to the discreetness of {ti — t2)/2, i.e. equal to a half of 
the width of one TDC channel Alps/ 2. 

To cross-check the results, another method was employed as well. The 
center of the peak gravity Mpeak and the mean square deviation apeak were 
directly calculated from the sample of collected events as 

Mpeak = (16) 



^peak 



\ ^H^Mpeak - X,)\ (17) 



where Nev denotes the total number of selected events, Xj is the time dif- 
ference derived from the TDC readouts chtdcu ^'^d chtdc2^ for each recorded 
event as 

Xi = r{chtdcu - chtdc2,)- (18) 

Here r denotes the TDC scale (i.e. the width of one TDC channel). In 
the data analysis it was fixed to r = A7ps/ch. In reality r is affected by the 
TDC differential non-linearity. For the LeCroy2228B TDC unit used in these 
measurements, it was measured and found varying from 45 to 49ps/ch. This 
variation leads to systematic uncertainty in each measurement. The way to 
estimate and to reduce it is explained in the Section [HI 
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Figure 10: Example of (ti — t2)/2 spectrum. Gaussian fit of tlie peak is one 
of the ways to extract the TOF resolution. 
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The width of the peak is in fact the combination of the TOF and electronic 
resolutions. 

'^peak = \J<^TOF + '^electronics i^^) 

The electronic resolution (Jeiectronics mostly originates from the TDCs res- 
olution. It was estimated from the spectrum of y in which the same PM 
signal generates the common START for the TDC and the delayed STOP 
(Fig. [3]). In fact, 99% of events in this spectrum were located in a single chan- 
nel. The resulting value is (Jeiectronics ~ 6.9 ps. The corresponding corrections 
were implemented following Eq. [191 

Both methods were verified in simulations. The values ti and t2 were 
simulated event-by-event using random generators. The generator for ti im- 
itated the measured ti spectrum while the t2 generator was Gaussian. The 
width of the Gaussian distribution is equal to 2aTOF- 

The generated ti and t2 values were then digitized in accordance with the 
discrete scale of LeCroy2228B TDC and recorded in the same format as the 
experimental data files. The experimental and simulated data were analyzed 
using the same codes. The simulated TOF resolution was reconstructed with 
the accuracy of 0.1 ps in the case of direct calculation and 0.3 ps using the 
Gaussian fit. 

One well-known problem of timing measurements is the effect of time 
walks of constant-fraction discriminators (CFD). CFDs generate additional 
time shifts (time walks) which depend on pulse heights. To extract the TOF 
resolution, the events were always selected using the criterion 

L, < ^ {Ai - pedi){A2 - ped2) < (20) 

such that the pulse heights Ai and A2 of the selected events are nearly 
constant (for example, events from the beam peak). Therefore the influence 
of time walks on was minimized. 

8 TOF resolution with fine-mesh R7761-70 
photomultiphers in comparison with R2083 
PMs 

In this section we report a comparative measurement of the TOF resolutions 
of two similar counters. One counter was equipped with fine-mesh R7761- 
70 photomultiphers, and the other with R2083 PMs. Both counters were 
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made of 2 x 3 x 50cm^ Bicron408 scintillation bars. The pliotomultipliers 
were directly attached to the bar butts. The diameter of the photocathodes 
of R7761-70 photomultipliers, of 27 mm, covered only ~ 80% of the butt 
surfaces of the ends of the scintillator bar. Therefore the light collection in 
this counter was lower by factor 0.8 than in the case of R2083 PMs. 

The beam irradiated the counters at their centers in the direction per- 
pendicular to the counter axis. The TOF resolution of both counters was 
measured at the same conditions: after the beam adjustment the counters 
were attached one-by-one to the collimation system. The replacement of the 
counters took few minutes. During that short break the beam was off. How- 
ever, the beam parameters were kept fixed. The beam spectra and the count 
rates were similar in both measurements (Fig. [7]). To avoid the effect of the 
TDC differential non-linearity, the TOF resolution of each counter was mea- 
sured successively six times. The scheme of the first measurements is shown 
in Fig. [3l In the second measurement the signal of that PM which does not 
trigger the acquisition, was additionally delayed for 1 ns. In the third mea- 
surement the delay was 2 ns. Then the PMs were switched and the series of 
the measurements with additional 0, 1, and 2 ns delays was repeated. 

The data were analyzed in an identical way: the same codes was running 
over the data files corresponding to 12 (6 x 2) measurements. Only events 
from the beam-peak area 31.5MeV < L < 32.5MeV were selected to extract 
the TOF resolutions. 

The results were obtained by using both the direct calculation and the 
fitting of [ti — t2)/2 spectra. They are summarized in Table [2] and Table [31 
The TOF resolution obtained with R7761-70 PMs was in addition scaled by 
factor VO.8 (last columns of Table [2], [3]) in order to account for 80% light 
collection due to the smaller diameter of the R7761-70 photocathodes. Such 
"effective" TOF resolution corresponds to the same number of photoelectrons 
produced at the R7761-70 and R2083 photocathodes. Further the corrected 
R7761-70 TOF resolution is used for the comparison with R2083 PMs. 

The statistical error of each measurement was ~ 0.1 ps. The systematic 
uncertainty mostly arises from the TDC differential non-linearity. It can be 
estimated as the deviation of six measured data points from their mean value 

1 ^ 



19 



Measurement PMtrigger Additional (Ttofu^oss ^tofhttsi-to crTOFcormrei- 

delay (ns) 



1 


1 





18.85 


23.84 


20.74 


2 


1 


1 


19.09 


19.84 


17.26 


3 


1 


2 


16.14 


19.70 


17.14 


4 


2 





16.30 


21.68 


18.86 


5 


2 


1 


16.43 


22.56 


19.64 


6 


2 


2 


17.60 


18.59 


16.17 



Table 2: TOF resolutions (ps) of R2083 and R7761-70 PMs directly calcu- 
lated from the samples of collected events 



Measurement PMtrigger Additional aroFR^oaa crTOFfl776i-7o crTOFcorRTrei- 

delay (ns) 



1 


1 





18.25 


22.67 


20.28 


2 


1 


1 


19.02 


19.67 


17.50 


3 


1 


2 


15.76 


19.19 


17.16 


4 


2 





16.72 


20.92 


18.70 


5 


2 


1 


15.97 


22.00 


19.68 


6 


2 


2 


16.88 


18.33 


16.40 



Table 3: TOF resolutions (ps) of R2083 and R7761-70 PMs obtained using 
Gaussian fits. 
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AcT' 



TOFsyst 



1 ^ 

\\ g JK^TOF, - OTOF^ean)^ (22) 



The estimates for the systematic error in one measurement is AaTOFR2083syst = 
1.21 ps and AaTOFRriei-iosyst = 1-58 ps for the direct calculation method, 
and AaTOFR2083syst = 1-17 ps and AaTOFR776i-70syst = 1-38 for the Gaussian 
fit. The error for the average of 6 measurements should be scaled by The 
averages for the R2083 resolution and for the corrected R7761-70 resolution 
are 

crTOFR2083 = 17.4 ± 0.49 crTOFR77ei-70 = 18.3 ± 0.64 (23) 
for the direct calculations and 

O-TOFR2083 = 17.1 ± 0.48 (TTOFR7761-70 = 18.3 ± 0.56 (24) 

The ratio of the R7761-70 and R2083 resolutions is 

O"TOF_R7761-70 



O'TOFR2083 

and 

CTOFR7761-70 



1.05 ±0.066 (25) 



1.07 ±0.062 (26) 

O'tOFR2083 

for the direct calculation and for the Gaussian fits respectively. Further the 
average of both methods 

aTOFR77ei-70 ^ ^ ^ q q^^ ^27) 

O'TOFR2083 

is used. This ratio means that, if the number of photoelectrons is the same, 
the TOF resolution with R7761-70 photomultipliers would be ~ 6.6 ± 6% 
worse than in the case of R2083 PMs. 

The result proves the good timing performance of R7761-70 PMs. More- 
over, this ratio was obtained at low HV ~ 1300\^ of R7761-70 PMs. It is 
well known that PM timing properties becomes better at higher HV (see, 
for example, Ref. [6J). In addition, as it will be discussed in the section [TT| 
o"TOF_R776i-70 might be better in magnetic field. 
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9 Dependence of TOF resolution on light out- 
put 

To retrieve the dependence of the TOF resolution on the hght output, the 
sample of collected events was divided into many sub-samples following the 
criterion 



U < V(^i - Pedi){A2 - ped2) < U+i (28) 

The width of the bins Lj+i — Lj was 1 MeV. For each bin the width of 
the peak a^eak was derived using both methods. 

The results collected with R7761-70 PMs in the October 2008 are shown 
in Fig. [m The data obtained by two methods are consistent. The minor 
deviation at the lower light output is explained by the increasing contami- 
nation of scattered protons. This background differently affects the results 
obtained by each method. 

The data points are well fitted by The results of the fit are 

aTOF{j^) = ; (29 j 



'L{MeV) 

with x'^ = 2.4 for the directly calculated data points and 



97.2±0.36stat±2systips) 
lL{MeV) 



for = 0.95 for the Gaussian fit. The systematic uncertainty of 2 ps orig- 
inates from the TDC differential non-linearity and from the accuracy in the 
determination of the MIP position. 

10 Dependence of TOF resolution on coordi- 
nate and track angle 

The TOF resolution is related to the PM timing resolutions (Tpmi and crpM2 
as ^ 

^TOF = -^\J<^'pMl + ^PM2 (31) 
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Figure 11: Dependence of the measured time-of-flight resolution on the light 
output. Circles are directly calculated from the sample of collected events. 
Open squares are the results obtained using the Gaussian fit. Error bars 

correspond to statistical errors only. Where not visible, the errors bars are 
smaller than the symbol size. The curve is the fit of experimental data. The 
difference between two fits corresponding to each data set is not seen in the 
plot. 
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Figure 12: Dependence of the TOF resolution on the coordinate along the 
counter axis. 



where (JpMi and (JpM2 are the effective resolutions in each PM channel. If a 
scintillation occurs in the middle of a counter, the numbers of light photons 
that reaches each PM, are equal. Therefore ctpmi = 0"pm2- If a scintillation 
is located near a counter end apMi 7^ 0"pm2- Due to the exponential light 
attenuation inside a counter the number of photons which reach each photo- 
multiplier is different. This effect may generate the dependence of the TOF 
resolution on the axis along the counter axis. 

To retrieve this dependence, the TOF resolution was measured at different 
coordinates along the counter axis. The results are shown in Fig. [T2i The 
TOF resolution is ~ 15% better near the counter ends. 

In addition, two measurements were performed with the beam directed 
at 60° relative the counter axis. Within the systematic uncertainty, the TOF 
resolution was the same as that obtained with the perpendicular beam. More 
accurate data for the aroF dependence on the coordinate and the angle will 
be obtained during next beam runs. 
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Dependence Gain VS MagneticField | 
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Figure 13: Dependence of R7761-70 gain on magnetic field. 



11 Operation of R7761-70 PMs in magnetic 
field 

In this section we present first results on the operation of R7761-70 pho- 
tomultipliers in magnetic field. Both photomultipliers were placed inside 
air-cooled solenoids. The solenoids were designed and manufactured by the 
KNU group. They comprise two parallel sections each wounded with 8 layers 
of 1-mm dia cooper wire. The 5 mm gap between the sections allows the air 
passage in order to improve the cooling efficiency. In future the number of 
sections will be increased to 4. 

Each section generates 100 Gauss of the magnetic field per 1 A of the 
current. Due to heating, the maximum current is limited to ~ 6 A. Accord- 
ingly, the upper limit of the magnetic field in the two-sectional solenoid is 
1200 Gauss. Significant increase of the generated field should be expected if 
the air cooling would be replaced with the water-cooling system. 

The measurements were carried out first at Gauss, then at 1100, 1000, 
900 ...100, and again at Gauss. All the data takings were done at the 
same conditions one-by-one. The field was varied by changing remotely the 
solenoid currents. 

In total, 4 series of data were taken. In the first series the direction of 
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Dependence sigma_TOF VS MagneticField | 
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Figure 14: Dependence of R7761-70 TOF resolution on magnetic field. 



the magnetic was chosen parallel to the PM axis. In the second series the 
signal of one PM was additionally delayed for 1 ns. Then the direction of the 
magnetic field was turned to 180 deg by switching the polarity of the power 
supply. After that the 3rd and 4th series were done in the similar way. 

The results are shown in Fig. [13] and Fig. [TH The gain of R7761-70 PMs 
was attributed to the position of the beam peak. The beam spectrum in this 
measurement is shown in the panel 6 of Fig. HI No any dependence of the 
gain (i.e. shift in the peak position) on the magnetic field was observed. 

Surprisingly, the TOF resolution becomes better if the magnetic field is 
oriented parallel to the PM axis. The effect reaches ~ 8% at 1100 Gauss. 
The antiparallel field does not affect the TOF resolution. This observation 
has to be checked in next beam runs. 

12 CTOF Estimates 

In real CTOF assembly most of fast particles will be detected at forward 
angles. Below we present the estimates of the expected TOF resolution for 
minimum- ionizing particles emitted from a target at 90° and 45°. 

A minimum-ionizing particles passing through a 3-cm thick counter per- 
pendicular to its axis, deposit ~ 6 MeV of energy (Fig. [H]). The corresponding 
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TOF resolutions are (yTOFB,2Qsz = 17.3w31.5/6 ~ 39.6 ps and (TTOFt 



"JJ7761-70 



18.3^31.5/6 ^ 41.7 ps for R2083 and R7761-70 PMs respectively. 

The scintillation bars in the CTOF assembly will be viewed through light 
guides. The long bent light guides in the CTOF design with ordinary R2083 
PMs (Fig. [1]) deliver ~ 30% of scintillation light [Ij. The estimate for the 
expected CTOF resolution with R2083 PMs would be 

aTo™83LG(90°) = X 39.6ps ^ 72.3ps (32) 

In the CTOF design with fine-mesh photomultipliers (Fig. [T]), the light 
guides will shorter and not bent. One may assume their light transfer effi- 
ciency around 50%. The corresponding estimate at 90° is 

aroFi?776i-70LG(90°) = X Al.lps ^ 59.0ps (33) 

It is essentially better than that number for R2083 PMs. 

Minimum- ionizing particles emitted at 45°, deposit \/2 times more en- 
ergy in the CTOF counters than ones emitted at 90°. The expected CTOF 
resolutions are 

C"rOF_R,2083LG(45°) = (34) 

1 

0.85 ^ , X 39.6ps ^ 43.5ps 

0'TOFm7&l-7QLG{^^°) = (35) 

0.85 , X Al.lps ^ 35.4ps 

\/2ya5 

where the factor 0.85 originates from the dependence of the TOF resolution 
on the coordinate. 

The simulated particle identification (PID) based on the quoted above 
estimates is shown in Fig. [151 The particle speed [3 is used as an identification 
parameter 

where Itof is a distance from a vertex point inside a target to a hit point on 
a counter surface. The simulations include four sources of uncertainties: 
i) the TOF resolution and its dependence of the hght output; 
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Figure 15: Simulated particle identification in the CTOF detector with fine- 
mesh R7761-70 photomultipliers at 90° (left) and 45° (right). 



ii) the coordinate resolution = cftof x 15cm/ns; 

iii) CTOF azimuthal granularity; 

iv) The reconstruction of the vertex point ±1 mm. 

At 90° kaons are discriminated from pions up to ~ 600 MeV/c. Protons 
are discriminated from kaons up to ~ 1000 MeV/c. In the most critical 
region near 45° PID is essentially better: protons are separated from kaons 
up to ~ 1400 MeV/c, and pions from kaons up to ~ 850 MeV/c. 



13 Summary and future R&Ds at KNU. 

The measured ratio of the effective TOF/timing resolutions of R7761-70 and 
R2083 PMs ''TOFR77^i-7o ^ ^ ^ Q proves the advantages of the CTOF 

O"TOFH2083 

design with fine-mesh photomultipliers. This design will be more simple, 
less expensive, and will provide better performance than the "conservative" 
design with ordinary R2083 PMs. 

More information on the properties of fine-mesh photomultipliers will be 
obtained in the next beam runs. The R&D program includes: 

- measurements of the relative TOF resolution of several counters equipped 
with R7761-70, R5924-70 and R2083 PMs; 

- operation of R7761-70 and R5924-70 photomultipliers in magnetic field up 
to 0.4 Tesla; 

- accurate measurements of the dependence of the TOF resolution on the 
coordinate and track angle. 
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At the next stage silicon photomultipliers, micro-channel plates could 
be tested as well using the same technique. Another option might be the 
development of a time-of-flight system for the detection of neutrons. The new 
neutron beam line is currently under construction at the MC50 Cyclotron. 
It will offer a tool for testing prototype counters of the CLAS12 neutron 
detector [T7|. 
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